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Abstract 

We present a unified framework for the calculation of defect energies for those defects that can 
be represented as a superposition of isolated dislocations, and obtain both self and interaction 
energies of combinations of grain boundaries and cracks. We recover in special limits several well 
known quantities such as the energy of a low-angle tilt boundary, as well as other lesser known 
results, including boundary/boundary and crack/boundary interaction energies. This approach, in 
combination with simple dimensional analysis, permits the rapid calculation of defect energetics in 
the elastic limit. 
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I. INTRODUCTION 



The presence of extended defects in crystalline solids, such as grain boundaries and cracks, 
can dramatically affect their mechanical response. For example, in large-grained polycrys- 
talline materials, grain boundaries may sometimes impede dislocation motion with a con- 
comitant increase in yield strength,J*i while sharp cracks are stress concentrators that may 
initiate material fracture^. From a micromechanical point of view, the elastic fields as- 
sociated with grain boundaries or cracks lead to defect interactions that influence fracture 
behavior in highly-flawed systems and solute segregation^, thereby affecting strengthening 
or embrittling mechanisms. 

Calculations of extended defect interaction energies most often begin with the solution 
of the appropriate elastic boundary value problem, followed by the construction of the cor- 
responding energy density in terms of the stress (or strain) fields and elastic constants and 
the integration of this density over some volume of space. An alternative approach to the 
defect interaction problem follows from a micromechanical model of the defect and is based 
in the elastic Green function of the medium. Such an approach was pioneered by Mura 
and othera^, and will be applied here to defects that may be regarded, in some limit, as 
composed of elemental straight dislocations. As will be seen below, this formalism permits 
the straightforward calculation of interaction energies for different extended defects and, 
moreover, facilitates an intuitive, multipole-based analysis that reveals their asymptotic de- 
pendence on defect separation. 

We address in this paper two prototypical systems: grain boundaries and cracks, as these 
defects can often be modeled in terms of spatial distributions of line defects, the former 
typically by regular array(s) of straight dislocations each having a constant Burgers vector, 
and the latter by a continuous, localized distribution of dislocations with a Burgers vector 
density depends upon loading conditions. Thus, the eigenstrains corresponding to each 
defect are conveniently expressed in terms of a dislocation density tensor that embodies the 
relevant defect length scales and separations. For the purposes of illustration we consider the 
energetics of a low-angle tilt boundary, the interaction between tilt boundaries, the energetics 
of an array of cracks, and the interaction between an isolated crack and a grain boundary. 
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The application of this approach to other, related examples is then discussed. 

II. FORMALISM AND SELECTED ILLUSTRATIONS 

As indicated above, our calculations are facilitated by describing a defect (e.g., grain 
boundary, crack) in terms of a dislocation density tensor, . This solenoidal, second-rank 
tensor carries both local Burgers vector and line direction information, with a divergenceless 
condition that is a consequence of the topological constraint of line continuity. For example, 
the components of the density tensors corresponding to straight edge and screw dislocations 
at ri = r2 = with Burgers vectors 6, and line directions along are given, respectively, by 



Note that, in the first equation, it is assumed that the Burgers vector is aligned either along 
ri or r2, and so an arbitrary alignment in the plane may be constructed from appropriate 
linear combinations. 5{r) is the Dirac delta function. 

The elastic energy of the system E [p] can be written as a functional of p, as shown 
by Kosevich and others^ and, more recently, has been employed by Nelson and Toner— to 
investigate the effect of unbound dislocation motion on shear response in solids, by Rickman 
and Vinalsii to model the collective motion of dislocation ensembles, and by Rickman and 
LeSar^- to quantify the temperature dependent interaction of fluctuating dislocation lines. 
For our purposes, it is convenient to express the energy functional for an elastically isotropic 
medium with shear modulus p and Poisson ratio v as the Fourier integral 



where the integration is over reciprocal space (tilde denoting a Fourier transform), the kernel 
(without defect core energy contributions) is given by 








QikQjl + CilCkj + 



1 - u 



CijCki 



(3) 
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and Q and C are longitudinal and transverse projection operators, defined by 



Qij 


= Sij - 






Cij 


— ^ijl — 

q 




(4) 



respectivel}*iS, where 5ij is the Kronecker delta, and eijk is the Levi-Civita tensor. We 
note that it is possible to include defect core energies in a somewhat ad hoc fashion by 
augmenting this kernel by a term quadratic in the dislocation density^'^ A somewhat more 
realistic description of the core could be given by an energy term that depends on atomic 
coordinates, although we will not consider such core models as they fall outside of the 
mesoscopic description outlined here. As illustrated below, the distance dependence of the 
defect energy can be determined in a straightforward manner from Eq. ^ by employing a 
superposition of prototypical densities. 

A. Edge Dislocation Arrays - Low- Angle Tilt Grain Boundaries 

Consider first a linear array of N edge dislocations, each separated from its nearest neigh- 
bor by I and having a Burgers vector aligned along either the ri or r2 axes (J = 1 or 2, 
respectively), as shown in Fig. As is well known, this system, with J = 1 and in the 
limit — s> oo, is a model of a low- angle tilt grain boundary where, in the limit of small grain 
misorientation angle 0, £ ~ h/&^. From Eq. (j21) the dislocation density is 

pij{r) = b 5i3 5jj 5{ri) ^ 5(r2 - ni), (5) 

n 

where the summation is over all dislocations in the array. For convenience, we will work 
with the corresponding Fourier transform 

pij{q) = 27ib 5{qs) ^ exp {iq2ni) Si3 5jj. (6) 

n 

Before calculating the corresponding boundary energy, we note that singularities inherent 
in this continuum model that lacks a short distance cut-off lead to a divergent self-energy. 
The origin of the divergence is the fact that the energy density of a single edge dislocation 
diverges like 1/r, where r is the radial distance to the dislocation line. In reality, each 
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FIG. 1: a.) A low-angle, tilt grain boundary modeled as an array of edge dislocations. Each 
dislocation has a Burgers vector b and is separated from its nearest neighbor by £ ~ 6/0 for small 
misorientation angle ©. b.) Two such grain boundaries separated by a distance L. 

dislocation has an atomic core region that is not amenable to a continuum description. We 
introduce a small cutoff parameter a representing the core radius by replacing the real-space 
delta function by the broadened delta representatiooi^ 

^{r, - nl) = - ^—^ -. (7) 
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The corresponding Fourier transform can be obtaining by using the shifting property of 
transforms to obtain 

A{q2) = exp{iq2n£) exp{-a\q2\). (8) 
Therefore the elastic energy of the array of edge dislocations is 

Ej = 7T7^ / \ K,j,j{q) p,M P,.j{-q), (9) 



2(27r)3 J ^ g2 
where 

K^JzAq) = Q33QJJ + \^CsjCsj (10) 
and J = 1 or 2 (no summation) denotes the orientation of the Burgers vectors. More 
specifically, noting that q'^ = qf + q2, 

1-uJ qf + ql 

and the energy per unit length in the direction is given by 

ej = TTT^ r / / A Siq2) exp (-2a|g2|), (12) 



KsjsAql = ( ^ ) (11) 



(27r)2 (1-^)7 J {ql + qlf 
where a structure factor S{q2) has been defined as S{q2) = | Xln^^P (^9'2^^)P 
N{27T/i) Em^te - 27rm/£)i^ Explicitly we have 

N^ib^ / . , /"j (27rm/£)2 



ei = - — -— > exp (— 47ra|m|/£) / dqi „ , 

]\T 1,2 m=oo „ 2 

62 = , , ^ y exp(-47ra|m|/£) dq^ (13) 

Upon evaluating the first integral in Eq. (fT^ . one finds 

Nfib'^ exp (— 47ra|m|/£) 

= 2(2vr) (1 - u) H ' ^ ^ 



and therefore 



^1 = ~7TTn A ^ ~ -47ra/£ ^ In 

[Att) (1 — u) {AiTj (1 — u) 



(15) 



47ra 

where the latter approximation holds when a/i is small. If this array represents a planar 
defect (i.e., grain boundary) then the corresponding energy per unit boundary area is 



ei fib"^ 
^ Ni (47r£) (1 - z/) 
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Ana 



(16) 



This result has the same functional form as the classic Read-Shockley energy for a low-angle 
tilt boundary, though the constant in the argument of the logarithm depends on the details 
of the short-range (core) cutoff employed here-^. In this regard, we note that a similar 
approach to this problem for two-dimensional systems has been outlined elsewhere^. 

Finally we examine 62 for the extended defect for which b \\ r2. From Eq. (jlHj) it is clear 
that 62 — > 00 owing to the divergence associated with the m = mode. This result is expected 
since some components of the long-ranged stress field (namely (722) tend to a constant as 
X —* 00, and thus the volume integral over the energy density diverges in an infinite system. 
Given this behavior, it is evident that a low-angle boundary with this geometry cannot have 
dislocations with Burgers vector components in the boundary plane. 



B. Dislocation Array Interactions 



Consider next two low-angle, symmetric, tilt boundaries separated by a distance L, as 
shown in Fig. ^p. The interaction energy for these boundaries follows from the corresponding 
dislocation density tensor 



5{n) ^5(r2-n£) + 5{ri - L) ^5( 



r2 — mt 



(17) 



that can be transformed to yield 



exp {iq2ntj + exp (iqiL) ^^exp(ig. 



'2mt 



(18) 



The interaction energy for separated dislocation distributions is obtained from the cross 
terms (i.e., no self-energies) in Eq. 0. Following the procedure outlined above, one first 
calculates the interaction energy per unit length of dislocation line 



1,= — 00 v^i 



Upon evaluating this integral one finds that 



cos(giL) exp (— 47ra|m|/£) (19) 



Npb'^ 



2{\-v. 



m=l 



1 



2Ti\m\ 



exp (—a \m\ 



(20) 



where a = (Ana/i) + (27c L/i). The resulting summations can also be performed to yield the 
grain-boundary interaction energy per unit area 



lint 



/27rL 



In [1 — exp {—a)] 



(21) 



Ni 27v{l-u)£ [\ i J Vexp(a)-1. 
It is of interest to examine 'jint for two limiting cases. First, for large boundary separations 
such that L/£ » 1 and a/i << 1, Cmt oc (L/i) exp {—2tiL/€). This exponential decay in the 
interaction energy follows from the rapid (exponential) decay of the stress fields associated 
with the individual arrays. By contrast, for L/l ~ 1, the logarithmic dependence of Cint 
on L/^ results from individual dislocations in either array "seeing" each other. Finally, we 
note that an alternative route to Eq. ()21|) follows from a calculation of the Peach-Koehler 
force acting on a dislocation owing to a distant array and a subsequent spatial integration 
to obtain the energy.^^^ 



C. Crack Array 

Another application of the formalism presented above is the interaction energy between 
cracks. We first consider a single crack of length 2c oriented along the ri axis, and model it as 
a continuous dislocation distribution, with corresponding dislocation function i?(ri).— The 
Burgers vectors of the dislocations that comprise the crack model embody the local crack 
opening that results from a given loading. Hence, by choosing an appropriate dislocation 
distribution, one can represent loading in various modes. Having represented the crack with 
a dislocation distribution, the stress field associated with the crack is given in terms of the 
stress field of a single edge dislocation a^^ by the convolution integral 

/ + 00 
dr[B{r[)al-^{n-r[,r2) (22) 
-oo 

As a specific example, consider a Mode II crack. The shear loading associated with this mode 
can be represented by a distribution of edge dislocations with Burgers vectors oriented along 
ri. Thus, one can write Biji) = J dr2 Pziiji,'r2) for an appropriate dislocation density psi 
that can be determined from the requirement that the crack faces must be traction-free.-^ 

It is of interest here to consider a crack interacting with distant objects. In particular, 
given an observation point r such that |r|/2c >> 1, it is permissible to regard the crack fields 
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FIG. 2: An array of cracks, each separated from its nearest neighbor by a distance i along the 
r2 axis. It is assumed that the system is subjected to a constant shear stress r. Also shown is a 
schematic of the dislocation distribution corresponding to each crack. 

as produced by the lowest-order multipole moments of B{ri). In practice, the "monopole" 
and "dipole" moments are often satisfactory in this context and one obtains 

/ + 00 
dn Bin) = btot, (23) 
-oo 

and, 

D = j dn n Bin) = — — — - J dn \Jc^ - rf ai2in,r2 = 0), (24) 

where btot is the total Burgers vector associated with the crack dislocations and 0"i2(ri, r2 = 0) 
is the shear stress loading the crack . 

We now derive the interaction energy associated with various crack arrays. Consider, for 
example, a linear array of Mode II cracks, each crack closed at both ends (i.e., btot = 0) 
and modeled as a dislocation dipole in the r2 direction that is separated from its nearest 
neighbor by i. As is evident from Eq. (j23), and by contrast with the grain-boundary models 
discussed above, the strength of elemental defects (i.e., their dipole moment) is not fixed, as it 
depends on the local stress field. The geometry for this loading, as well as the corresponding 
dislocation orientations, is shown in Fig. |21 In the short crack limit, such that each crack 
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may be regarded as a point dipole with moment Di (i.e., the hmit c ^ and h ^ oo with 
2ch — > Di), the corresponding dislocation density tensor is 



P31 



{r) = DAs S,i S'in) J2 '^(^2 - ni), (25) 



where the prime denotes differentiation with respect to the argument, and the subscript i 
indicates that the dipole moment depends on craclc separation. As before, it is convenient 
to work in reciprocal space where one finds that 

Psiio) = 27iiDe qi 5{qs) ^exp {iq2n£) Sis 5ji (26) 

n 

By analogy with the development given above, the energy of this system is given by 

NfiDj ^ f, 2 (27rm/£)2 ^ ^ , , 

(27r) (1 - v)l J (g2 + (2vrm/£)2) 

Upon performing the required integral and summation, one finds that 

_27r NfiDj exp(47ra/£) 

= (1 - u)i2 [exp (4vra/^) - if ^^^^ 

It is of interest to examine Cc in the limit of large crack separations i where the dipole approx- 
imation works best. One finds that [NfiDj /8n{l - u)] - An^/Sf + 0(1/ 1)], the 
first term arising from the self-energy of a dislocation dipole, and the second arising from 
dipole-dipole interactions. 

An approximation to Di for a particular array spacing i may be obtained by requiring 
that the shear stress traction produced by neighboring cracks on the face of a given crack be 
compensated by a uniform stress that is the same for every crack in the array. The strategy 
for such calculations is given elsewhere^^, and the details for this defect geometry are given 
in Appendix 1X1 



D. Crack-Grain Boundary Interaction 

As a final illustration of extended defect energy calculations, consider the interaction of 
a low-angle grain boundary separated from an isolated crack by a distance L, as shown 



10 



2c 

< 



i 



FIG. 3: An isolated crack interacting with a low-angle tilt boundary, the latter modeled as an array 
of edge dislocations. For simplicity, the crack is located at a distance L from the boundary and at 
r2 = 0. For this geometry the crack is subjected to shear stresses owing to the grain boundary. 

schematically in Fig. IHl For this geometry the normal stress on the crack faces owing to the 
boundary is o"22(?"i, ^2 = 0) = 0, while the shear stresses au^ri, r2 = 0) 7^ 0, and so the crack 
is loaded in Mode II. From the grain-boundary (Eq. ©) and crack (Eq. dislocation 
densities one can then construct the interaction energy per unit dislocation length, namely 



''9b = Y] \2 (l^ — V / / , n \,2 ^ 9i (^9i^) COS {(i2nt) exp (-a|g2|) (29) 
n ^^'^> J J {q{ + g|) 



Evaluating the integrals and the summation, we find that 

a + L 



_ fibPL Lit'' 2 
""^^^ " (27r)2 (1 - u) P 



TT 



(30) 



L 

In the limit of i/L << 1, Ccgb — {L/i"^) exp [— 27r (a + L) /L], owing to the exponential decay 
of the stress fields associated with the grain boundary while, for £/L ~ 1 (and a/L << 1) 
Cint — '^/L, consistent with a monopole-dipole interaction. We note that this result can 
be obtained somewhat more readily by returning to the calculation of the grain boundary 
interaction energy and considering a single dislocation within a grain boundary array. Since 
a crack is modeled here as a dislocation dipole one replaces, in effect, an edge dislocation 
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density with its spatial derivative and therefore the functional form of the crack/boundary 
interaction energy is given by L (deint/dL) (see Eq. fET|)). 

The crack dipole moment can again be determined from the shear loading conditions 
via Eq. ()24j) . For the crack geometry shown in Fig. IHl this leads to 

-1 



-2bn 



dri \ — r\ (ri + V) 



cosh 



27r (n + L) 



1 



(31) 



where the shear stress is that for an array of edge dislocation.^^ This integral can be evalu- 
ated numerically by a Gauss-Chebyshev integration technique based on type II Chebyshev 
polynomials.— 



III. CONCLUSIONS 



A unified framework for the calculation of defect energies has been presented here and, for 
the purposes of illustration, applied to several different systems containing grain boundaries 
and/or cracks. In each case, a set of elemental dislocations comprise the extended defect 
under consideration, and the corresponding dislocation density is either independent of other 
defect positions (e.g., grain boundary) or calculable from the stresses imposed by these other 
defects (e.g., crack). As indicated above, this approach is especially useful in providing an 
intuitive understanding of interactions based on idealized defect models. 

The foregoing development can, of course, be applied to other systems for which a 
dislocation-based model is appropriate. For example, one can also obtain the energetics 
of a roughened tilt boundary^SiSl by considering perturbations of the dislocation density re- 
sulting from sinusoidal variations in the position of constituent dislocation lines. If these 
variations in boundary morphology are temperature induced, this analysis can be used to 
compute the statistical weights of perturbed boundary configurations to the free energy of 
the system and, consequently, to deduce a thermodynamic roughening temperature. In ad- 
dition, the description of more complex, asymmetrical boundary structures, consisting of 
two or three sets of edge dislocations, or twist boundaries consisting of perpendicular sets of 
screw dislocations is also possible by a generalization of the structure factor (see below Eq. 
(I12|) ) to include a form factor that reflects the positions of the basis dislocations within a 
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unit cell that generates the boundary. These and other calculations will be the subject of a 
future publication. 

Finally, we mention that, for the sake of completeness. Appendix El contains expressions 
for the stress tensor in the same Fourier representation used here for the energy, along with 
some simple examples. In a future publication, we will apply a similar formalism as used for 
the energy to derive equations for the stress arising from extended defects. 
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The Green function of the biharmonic operator can be obtained by noting that, in infinite 
space, V^|r| = V^V^|f| = —6{r). This result follows from the Green function relation for the 



three-dimensional Laplace operator V^|r[ 



1 

'FT- 



APPENDIX A: CRACK ARRAYS 



An approximation to the stress on a crack needed to balance those owing to others in the 
linear array shown in Fig. |21 can be obtained via the formalism of Dyskin and Miihlhaus.-^ 
We follow their approach below, except that the required stress fields are obtained by a 
multipole expansion rather than by using the Muskhelishvili complex potentials.— First, 
assuming that the crack separation i is large, one can regard each crack as a dislocation 
dipole. For this geometry the relevant stresses associated with each crack, in the limit of 
small crack width 2c, are given in terms of the derivatives of the derivative of the stresses 
(i.e., a point dipole approximation) for individual dislocations by 



0-12 (ri,r2j 
and 

0-22 (n,^2) 



libc 



d 



TT (1 — z/) dri 



jibe d 



ri {rl — rl) 



2n(l - u) 



3rl — rl 



Arl {rl — r|) 

2^ 



TT (1 — z/) dri 



r2 {rj 



.2\2 



fiD rir2 



2 {rl - rl) 



.{rl + rlf {rl + rlY_ 



2^3 



(Al) 



(A2) 



where D is the dipole moment. Given the symmetry of the crack array, it is clear that the 
stresses (T22 will cancel upon summation over all cracks. 

Next, we invoke a dipole asymptotics approximation^^ in which the i-th crack is subjected 
to a loading shear stress r and an additional uniform load e^, the latter equal to the stresses 
generated by the other cracks at the center of the z-th crack. The corresponding dipole 
moment is Di = — ttc^ (1 — z/) (r + e^) /fi. The corrective stress associated with the j-th 
crack is then given by 



— y 



1 



(r + e,) . 



(A3) 



Finally, taking each crack to be identical so that = e and using the Riemann zeta function 
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relation l/in? = ({2) = 7r^/6,— one finally obtains 



e^r '-/^''f' . (A4) 



The dipole moment can now be expressed in terms of r and e. 



APPENDIX B: STRESS TENSOR 



The stress tensor can also be given in terms of the dislocation density tensor. As discussed 
by Kosevicb^, the stress field at point r 



1 — z/ \dxidxk 



with 



and 



on 



Vik = - 



dpki dpii_ 



In Fourier space one can write Eq. ()B2|) as 



(Bl) 



:b21 



(B3) 



l<'^k{q) = G{q) fj^kiq), (B4) 
with G{q) = 1/g^, the Green function of the biharmonic operator--. Furthermore, 

Vikiq) = - [(^ipiqpPiiiq) + ekpiqpPiiiq)] = -— [Cupkiiq) + Ckipu^q)] ■ (B5) 

Given these transforms, the Fourier transform of the stress tensor can be written, after some 
algebra, as 



(^ik{q) = ipq^G{q) (Ciipki + Ckipii - 

V 1 - 



QikCalP. 



al 



(B6) 



To illustrate the use of Eq. ()B6|) . consider the case of a single straight screw dislocation, 
aligned along and having Burgers vector b. The corresponding dislocation density is 
Pki = 2Txh 6k3Si3 Siqs). 

Upon substituting this density into Eq. ()B6|) one finds that 



aik{q) = 2'Kipq^ G{q) b ^(ga) {Ci^Sks + C^s^is) • 



(B7) 
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Therefore, one immediately sees that an — 022 — eras = 0, and that a\2 — 0. The non- 
vanishing stress component 



ai3 = -i27r//6 ( % \ 5{qs). (B8) 



and so, in real space. 



-13 = (B9) 



A similar result can be derived for (723. 
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